Atomic Quantum Simulation of U (N) and SU (N) Non-Abelian Lattice Gauge Theories 
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Using ultracold alkaline-earth atoms in optical lattices, we construct a quantum simulator for 
U (N) and SU(N) lattice gauge theories with fermionic matter based on quantum link models. These 
systems share qualitative features with QCD, including chiral symmetry breaking and restoration 
at non-zero temperature or baryon density. Unlike classical simulations, a quantum simulator does 
not suffer from sign problems and can address the corresponding chiral dynamics in real time. 



Introduction. Non-Abelian gauge fields play a central 
role in the dynamics of the Standard Model of particle 
physics. In particular, the strong SU(3) gauge interac- 
tions between quarks and gluons in QCD give rise to 
the spontaneous breakdown of the chiral symmetry of 
the light quarks. Heavy-ion collisions produce a high- 
temperature quark-gluon plasma in which chiral symme- 
try is restored. The deep interior of neutron stars con- 
tains high-density nuclear matter or even quark matter, 
which may be a baryonic superfluid or a color super- 
conductor pQ. Unfortunately, due to severe sign prob- 
lems, the real-time evolution of heavy-ion collisions or 
the phase structure of dense QCD matter is inaccessi- 
ble to first principles classical simulation methods. In 
condensed matter physics strongly coupled gauge the- 
ories play a prominent role in strongly correlated sys- 
tems. In particular, the non-Abelian SU(2) variant of 
quantum spin liquids has long been debated as a pos- 
sible connection between the doped Mott insulator and 
the high-T c superconducting phase in cuprates jj,- The 
challenge of solving such problems motivates the devel- 
opment of quantum simulators for non-Abelian lattice 
gauge theories. Recently, quantum simulators have been 
constructed for Abelian U(l) gauge theories with [3] and 
without coupling to matter fields [H0. Here, we con- 
struct a quantum simulator of U(N) and SU(N) strongly 
coupled lattice gauge theories in (1 + 1), (2 + 1), and 
(3 + 1)D using ultracold alkaline-earth (AE) atoms in an 
optical lattice. On the one hand, our approach is based 
on quantum link models (QLMs) [MS]: which allow the 
exact embodiment of non-Abelian gauge interactions in 
ultracold matter. On the other hand, we utilize funda- 
mental symmetries of matter, such as the SU(2I + 1) 
invariance of interactions between fermionic AE isotopes 
such as 87 Sr or 173 Yb [SHU]- While still being far from 
a quantum simulator for full QCD, simpler model sys- 
tems share several qualitative features, including confine- 
ment, chiral symmetry breaking (^SB), and its restora- 
tion (%SR) [1 . They provide a unique environment to 
investigate important dynamical questions which are out 
of reach for classical simulation. 

The non-perturbative physics of non-Abelian gauge 
theories is traditionally addressed in the context of Wil- 




FIG. 1. [Color online] a) (upper panel) U(N) QLM in (1+1)D 
with quark helds tj/ x on lattice sites and gauge fields U^ x+1 
on links; (lower panel) hopping of AE atoms between quark 
and rishon sites of the same shading, b) Implementation of 
the QLM in rishon representation with fermionic atoms in 
(2 + 1)D. c) Encoding of the color degrees of freedom for 
N = 2 (t , 4) in Zeeman states of a fermionic AE atom with 
/ = 3/2. d) Lattice structure to avoid the interaction in 
fermionic matter sites using a species-dependent optical lat- 
tice (for an alternative method using site-dependent optical 
Feshbach resonances see the main text), e) Initial state loaded 
in the optical lattice with a staggered distribution of doubly 
occupied sites for a (7(2) QLM with M = 2. 

son's lattice gauge theory [15], in which the gluon field 
is represented by parallel transporter matrices residing 
on the links connecting neighboring lattice points of a 
4D space-time lattice. Since Wilson's classical link vari- 
ables take values in the continuous gauge group SU (N) , 
the corresponding Hilbert space is infinite-dimensional 
even for a single link. The elements of a quantum link 
matrix are non-commuting operators, similar to the com- 
ponents of a quantum spin. As a result, QLMs have a 
finite-dimensional Hilbert space, and therefore provide 
an attractive framework for the construction of quantum 
simulators for dynamical Abelian and non-Abelian gauge 
theories. In the continuum limit of QLMs, which is nat- 
urally realized via dimensional reduction, one recovers 
QCD with chiral quarks as domain wall fermions [201 [H] . 
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A pedagogical introduction to QLMs is contained in the 
supplemental material [22] . 

Quantum Link Models. The hopping of electrons be- 
tween lattice sites x and y in an external magnetic back- 
ground field B — V x A is described by tp x u xy ip y , where 
u xy — exp(i f v dl ■ A) £ U(l) is the phase picked up in 
this process [23]. In particle physics, gauge fields ap- 
pear as dynamical quantum degrees of freedom, not just 
as classical background fields. Here we consider U(N) 
and SU(N) lattice gauge theories without approaching 
the continuum limit, using so-called staggered fermions, 
which are represented by creation and annihilation oper- 
ators ip z J and ip x , that obey standard anti-commutation 
relations. Here i £ {1,2, ... ,N} represents the non- 
Abelian color index of a quark. The fundamental gauge 
degrees of freedom representing the gluon field are N x N 
matrices U xy (with elements U x ^ y ) associated with the link 
between nearest-neighbor points x and y (c.f. FigjlJ,). 
The hopping of a quark, which exchanges color with the 
gluon field, is then described by ^ x U xy ^y = ip£U* 3 y ipL 
This term is invariant against gauge transformations, 

n ^ x = n x ip x , n v4 = ^l n u xy = n x u xy nl with 

fl x £ U(N). The SU(N) gauge transformations and the 
additional U(l) gauge transformation contained in U(N) 
are generated by 

^ x ^ x ^ij^x ^ ; (^^x,x+k ^x~k,xj ' 
k 

g* = v4V*-E(^,* + *-^-m)> « 

k 

where A: is a unit-vector in the fc-direction, A a (a € 
{1, 2, . . . , A 2 — 1}) are the Gell-Mann matrices, and f a b c 
are the SU(N) structure constants, such that [G x , G y ] = 
2iS xy f a bcG x . The operators L xy and R x represent 
SU(N) electric fluxes associated with the left and right 
end of a link (xy), while E xy represents the Abelian 
U(l) electric flux. Physical states |*) obey the SU(N) 
Gauss law G X YS>) = 0, while in a U(N) gauge theory also 
G x \9) = 0. The operators U, L a , R a , and E associated 
with the same link obey 

[L a ,L b ] = 2if abc L c , [R a ,R b ] = 2tf abc R c , 

[L a ,R b ] = [E,L a ] = [E,R a ] =0, 

[L a , U] = -X a U, [R a , U] = UX a , [E, U] = U, (2) 

while operators associated with different links commute. 

In Wilson's lattice gauge theory, U is an element of 
the gauge group. In case of a U(N) gauge theory, 
dett/ = exp(iip) e [/(I) represents a C/(l) link vari- 
able, canonically conjugate to the electric flux operator 
E = -id v . In an SU(N) gauge theory U € SU(N) 
and L a , R a take appropriate derivatives with respect to 
the matrix elements U %: > . The resulting Hilbert space 
per link is then unavoidably infinite-dimensional. In or- 
der to represent the commutation relations of the gauge 



algebra of Eq.([2]) in a finite-dimensional Hilbert space, 
QLMs give up the commutativity of the matrix ele- 
ments U %0 . The real and imaginary parts of the ma- 
trix elements of the N x N quantum link matrix 
are represented by 2N 2 Hermitean operators. Together 
with the electric flux operators L a , R a , and E these are 
2A 2 + 2(A 2 - 1) + 1 = (2A) 2 - 1 generators which form 
the embedding algebra SU(2N). While U(l) quantum 
links can be represented by quantum spins embedded in 
an SU{2) algebra, U(N) or SU(N) QLMs can be realized 
with different representations of SU (2N) . A very useful 
representation is based on fermionic rishon constituents 

m 

u ij = cj^l, n = civ: + 4^. (3) 

The rishon creation and annihilation operators, c± and 
c±, are associated with the left and right ends of a link 
(c.f. FigjlJ,) and obey standard anti-commutation rela- 
tions. Our construction of a quantum simulator for U(l) 
gauge theories used Schwinger bosons to represent quan- 
tum links [3]. Here it is natural to replace Schwinger 
bosons by rishon fermions. J\f counts the number of ris- 
hons on a link. 

The Hamiltonian of a (d+ 1)D U(N) QLM with stag- 
gered fermions takes the form 

# = {s*v^U%% + h.c.) + m^ Sx vM 

(xy) x 

= -t E (s xv Q x+k Q y - k + h.c J +m^2 s xM x , (4) 

(xy) x 

where s x = (-1)^+-+^ and s xy = , 
with y = x + k. We have also introduced the U(N) 
gauge invariant "meson" and "constituent quark" oper- 
ators M x = ^ip x and Q x> ±k = c x ±k ip x . Together with 

the "glueball" operators $ x ,±k,±l = ej ±k c x ±i> they form 
a site-based U(2d + 1) algebra. The rishon number is 
conserved locally on each link. The U(N) model has 
no baryons, since the U(l) baryon number symmetry is 
gauged. In order to obtain charge conjugation invari- 
ance C and to reduce the gauge symmetry to SU(N), 
one must work with J\f xy = N rishons per link. Adding 
the term 7 Ylixy) (detU X p + h.c.) to the Hamiltonian, ex- 
plicitly breaks the U(N) gauge symmetry down to a lo- 
cal SU(N) and a global U(l) baryon number symme- 
try generated by B = J^ x yPx^x ~ t)- The symme- 
tries of various model systems are summarized in Table 
1. All models have a Z(2) chiral symmetry, which is 
spontaneously broken at a critical temperature T c , and 
may get restored at non-zero baryon density ns- It 
would be natural to add electric and magnetic field en- 
ergy terms ^- J2( xy ) (L xy L xy + R% y Rxy)i ^~ ^2( xy ) ^xyi 
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dimension 


group 


\r 
J\ 


o 


flavor 


baryon 


phenomena 


(1 + 1)D 


[7(2) 


1 


no 


no 


no 


XSB, T c = 


(2 + 1)D 


U(2) 


2 


yes 


Z(2) 


no 


X SB, T c > 


(2 + l)D 


SU(2) 


2 


yes 


Z(2) 


C7(l) 
boson 


XSB,T C > 
xSR,n B > 


(3 + l)D 


SU(3) 


3 


yes 


Z(2) 2 


C/(l) 
fermion 


X SB,T c >0 
xSR,n s > 



TABLE I. Symmetries and phenomena in some QLMs. 



and 4^12(wxyz) (U wx U xy U yz U zw + h.c), where (wxyz) 
denotes an elementary plaquette. At strong coupling 
these terms are inessential for qualitative features of the 
dynamics at finite temperature or baryon density, and 
are thus not yet included in our implementation. 

Atomic quantum simulation of U(N) QLMs. An il- 
lustration of the QLM and its rishon representation for 
(1+1)D and (2+l)D is provided in Figjl] Quark fields r/£ 
reside on the lattice sites x, while the rishons c z x ±k are 
on "link-sites" (a;, ±/s) at the left (right) end of the links 
exiting (entering) the point x (c.f. Figfl^i lower panel). 
The key step in our physical implementation is to inter- 
pret the lattice with quark and rishon sites in FigsjT^,,b 
as a physical optical lattice for fermionic atoms. Hence, 
an atom on site x of the optical lattice represents a quark 
tpl, while hopping of this atom to a link-site (x, ±fc) con- 
verts it to a rishon cj. ±k . The color index i is encoded 
in internal atomic states. 

The basic building blocks in our atomic setup are the 
tunnel-coupled triple wells in (1+1)D (FigjTa) or the 
cross-shaped vertices in (2+l)D (Fig[]J>). The corre- 
sponding hopping dynamics of the atoms is described by 
the Hamiltonian h x< k = t(s xy Q X}+ k+Q x ,-k+h.c.). Phys- 
ically, the overlap of the wannier wave functions can be 
used to implement the usual tunneling (24]. In case dif- 
ferent phases are needed to simulate staggered fermions 
in the lattice, Raman assisted tunneling 23J or shaken 
optical lattices [551 HI] can be applied. In order to ob- 
tain the desired quark-rishon dynamics, we introduce the 
microscopic atomic Hamiltonian 

H = U^{Af xy -n) 2 + Y^ h x,k + mY^s x M x . (5) 

( X y) x,k x 

The first term enforces the constraint of Af xy = n rishons 
per link, with U ^> t. In a physical setup this is imple- 
mented as a strong repulsion between atoms occupying 
rishon-sites, indicated in Fig{T^, by the overlapping link- 
sites. The second term represents atomic hopping, while 
the last term realizes the staggered fermion mass with a 
superlattice. In second order perturbation theory in the 
tunnel-coupling, the above Hamiltonian induces the hop- 
ping term of Eq.Q with t = P/U. FigjSji. illustrates the 
matter-gauge interaction. We note that an additional 
term t J2 X ±k Q x ±kQx,±k is also generated. This is no 
problem, because this term is invariant under all rele- 



vant symmetries. It is also possible to add a 4-fermion 

termFE^r 

With the Hamiltonian of Eq.d5j| we have reduced the 

realization of U(N) QLMs to a lattice dynamics of inter- 
acting fermions. This is enabled by the factorization of 
the quantum link variables into rishons. We emphasize 
that the building blocks in H are gauge invariant "me- 
son" and "constituent quark" operators, which allows a 
gauge invariant implementation of the dynamics. This is 
in contrast to previous work, where Gauss' law was en- 
forced by an energy constraint in the microscopic dynam- 
ics. The essential symmetries of H to be respected by the 
implementation are: (i) the color-independent hopping 
of fermions and rishons, and (ii) the color-independent 
interaction between rishons to ensure the local particle 
number conservation in each link. Indeed these symme- 
tries are accurately respected in setups with AE atoms 

EH mug. 

For a given nuclear spin /, the electronic ground state 
1 S'o of fermionic AE atoms has 21+1 Zeeman levels 
mi = —I, . . . , +1. We encode the color degrees of free- 
dom for the even (odd) building blocks (triple-wells in 
(1+1)D and cross-shaped vertices in (2+l)D, represented 
by the light (dark) shading in Fig[T]) in the N lowest 
(highest) mi levels (c.f. Fig[Ij:). For example, to imple- 
ment a U (2) QLM, we choose positive nuclear spin states 
mi = 3/2,1/2 on the even and negative nuclear spin 
states mi = —3/2, —1/2 on the odd building blocks. The 
AE atoms have the unique property that their scattering 
is almost exactly SU(2I + l)-symmetric, i.e., all pairs of 
states have the same scattering length [HJ [18] . This 
guarantees the symmetry of the U term in Eq.([5]). The 
mj-dependent hopping illustrated in Figflji can be real- 
ized in optical lattices with an appropriate choice of laser 
frequencies and polarizations [13 [55] , or with optical po- 
tentials obtained by holographic techniques [29l |30] . Fi- 
nally, the repulsion U, which only affects the rishon- but 
not the quark-sites, can be realized with optical Fesh- 
bach resonances of AE atoms allowing spatially depen- 
dent on-site interactions J31 36 . An alternative setup 
uses m/-dependent optical lattices with overlapping sites 
for the interacting, and spatially separated sites for the 
non- interacting fermions (c.f. FigjTJi) . 

SU [N) lattice gauge theories. We now reduce the 
gauge symmetry from U(N) to SU(N) by activating the 
detU xy term. For definiteness, we investigate the N — 
H = 2 case, for which detU xy = 24 )+fc e^_ fe e^ +/ .c^_ fe 
This corresponds to two-particle tunneling between the 
overlapping rishon-sites (c.f. Figj2]D). In our AE imple- 
mentation this process can be realized by Raman transi- 
tions coupling the internal states, with an intensity and 
a large detuning S, such that single particle transitions 
are excluded (see [37]). The two-body process is allowed 
by the fact that we explicitly break the SU(2I + 1) sym- 
metry of the AE setup down to SU(2) x ^ + k x SU(2) y -k, 
resulting from different overlaps of the Wannier functions 
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b) 



X + 1 



4- 



•Zj+^+l,- 



det U x . 



— lt)U)-U)lt) 

FIG. 2. [Color online] Dynamical processes in U(2) QLMs 
with N = 2. a) Matter-gauge interaction as correlated hop- 
ping of quarks and rishons. Starting with a configuration of 
site-singlets, the matter-gauge interaction converts them into 
nearest-neighbor singlets, keeping the rishon number per link 
constant, b) The determinant term corresponds to two-body 
hopping of both rishons on the link. 



in the partially overlapping rishon-sites. This generates a 
small difference AJ7 of the repulsive interactions between 
rishons on the same and on different link-sites. The re- 
sulting coefficient of the AetU xy term, 7 = 2Q 2 AU/S 2 , 
although small, should be larger than the typical tem- 
perature scale in cold atoms experiments. 

Initial conditions and loading the optical lattice. In 
order to evolve from an initial state that obeys Gauss' 
law, we now discuss how to load the lattice with the 
gauge invariant state illustrated in Fig{T^. This state has 
local color-singlet pairs of atoms on alternating quark 
and rishon sites. It is an eigenstate of the Hamiltonian 

^strong =™^2 S ^xi ) i + 7 ( detJ7 ^ + h - C -) ' ( 6 ) 

x (xy) 

which is induced by H in the limit U — > 00. The 
detU xy term favors the state | tli 0) — \0,ti) where 
I tl,0> = <£+c£! + |0 > 0) and |0,U> = c^_cf_ |0,0>. The 
preparation of the initial state requires to adiabatically 
ramp up the optical lattice on an ultracold cloud of 
atoms which are internally in a 50% mixture of the states 
mi = 3/2, 1/2. This leads to a band insulator with two 
atoms of positive nuclear spin on the dark-shaded sites 
111 Fig[l] Then, an on-site Raman two-body process will 
generate the desired state of FigflJ: after a coherent trans- 
fer of the rishon population from the dark- to the light- 
shaded rishon-sites. 

Imperfections and quality of gauge invariance. The 
emergence of gauge symmetry from the microscopic 
model of Eq.([5]) and the effect of imperfections have been 
addressed on small system sizes by means of exact diag- 
onalization in the U(2) case with Af xy = 2. The U(l) 
symmetry generated by G x has been verified by check- 
ing the condition (N xy ) = 2 + An n k(C//t, m/t), where 
Aii„k has to be small in order to ensure gauge invariance. 
As in Ref. [3], we found that a relatively small U is suf- 
ficient to accurately realize gauge invariance [22 . The 
SU(2) gauge symmetry, on the other hand, is potentially 
vulnerable to imperfections in the engineering of the tun- 
neling rate. Fortunately, symmetric interactions emerge 




FIG. 3. [Color online] a) Size L dependence of the energy 
splitting between the lowest energy eigenstates of a U(2) QLM 
with m = and V = —6t. b) Real-time evolution of the order 
parameter profile (iI>iI>)x(t) for L = 12, mimicking the expan- 
sion of a hot quark-gluon plasma. Here, circles (thin line), 
diamonds (dashed line) and squares (dotted line) correspond 
to r/t — 0, 1, 10 respectively. 



naturally thanks to the fundamental SU (21 + 1) symme- 
try of AE atoms. Since the number of fermions per ba- 
sic building block (including sites of the same shading in 
FigJT]) is conserved in the experiment up to exponentially 
small long-range tunneling processes, G x automatically 
commutes with H. One still needs to check whether the 
expectation values of the operators (G x ) 2 , (G 2 ) 2 vanish 
in the ground state. We have verified this numerically 
on small systems with up to 12 fermions (4 quark- and 8 
rishon-sites) by introducing an imperfect color-dependent 
hopping rate t\ = (1 — z)t 2 . We observe that gauge in- 
variance is very well preserved even for z ~ 0.1 |22j . 

Exact diagonalization results. We have performed ex- 
act diagonalization studies of the (1 + 1)D U(2) model 
with Af = 1 rishon per link. Figure [3^, shows the splitting 
between two almost degenerate vacuum states, which de- 
creases exponentially with the system size L, thus in- 
dicating 1(2) xSB. Figure [3|d shows the real-time evo- 
lution of the chiral order parameter profile (ipip) x = 
Sxty^ipx ~ %)> starting from an initial chirally restored 
"fireball" embedded in the chirally broken vacuum. This 
mimics the expanding quark-gluon plasma generated in 
a heavy-ion collision. 

Conclusions. We have proposed implementation of a 
quantum simulator for non-Abelian U(N) and SU(N) 
gauge theories for staggered fermions with ultracold 
atoms. The proposal builds on the unique properties 
of quantum link models with rishons representing the 
gauge fields: this allows a formulation in terms of a 
Fermi-Hubbard model, which can be realized with multi- 
component alkaline earth atoms in optical lattices, and 
where atomic physics provides both the control fields and 
measurement tools for studying the equilibrium and non- 
equilibrium dynamics and spectroscopy. Extending such 
investigations towards QCD requires the incorporation 
of multi-component Dirac fermions with the appropriate 
chiral symmetries, and of additional link and plaquette 
terms for electric and magnetic field energies [3"5] . 
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SUPPLEMENTARY MATERIAL I: U(N) AND 
SU{N) QUANTUM LINK MODELS WITH 
STAGGERED FERMIONS 

Introduction 

Non-Abelian gauge fields are the key to understanding 
the weak and strong interactions between leptons and 
quarks, mediated by W- or Z-bosons and gluons, respec- 
tively. The strong interactions are described by Quantum 
Chromodynamics (QCD), the relativistic quantum field 
theory of quarks and gluons with an SU (3) non-Abelian 
gauge symmetry. Due to asymptotic freedom, at very 
high energies quarks and gluons interact only weakly. At 
low energies, on the other hand, they interact so strongly 
that they do not appear as physical states in the spec- 
trum. Instead quarks and gluons are permanently bound 
together and thus confined inside hadrons. Protons and 
neutrons mostly contain up and down quarks as well as 
gluons, while the heavier strange quarks are more sup- 
pressed. Since the mass difference between up and down 
quarks is small compared to the intrinsic energy scale of 
QCD, there is an approximate isospin flavor symmetry 
SU(2)p. Protons and neutrons form an isospin SU(2)p 
doublet. In addition, the number of quarks is conserved 
which gives rise to a U{1)b baryon number symmetry. 
There are different categories of hadrons. States with a 
non-vanishing baryon number, i.e. with a different num- 
ber of quarks and anti-quarks, are known as baryons. In 
a simple-minded quark model, protons and neutrons are 
baryons consisting of three quarks of different SU(3) col- 
ors, while mesons are quark-anti-quark pairs. The light- 
est mesons of QCD are the three pions tt + , it , and w~ , 
which form an isospin SU(2)p triplet. The pions are 
the pseudo-Nambu-Goldstone bosons of the approximate 
SU(2) L x SU(2) R chiral symmetry of two flavor QCD, 
which gets spontaneously broken down to the isospin 
symmetry SU (2) p by the strong gauge interactions. Fla- 
vor singlet states with baryon number zero may be dom- 
inated by gluons and are then denoted as glueballs. 

Since quarks and gluons are confined inside hadrons 
(i.e. baryons, mesons, or glueballs) their low-energy dy- 
namics cannot be studied in perturbation theory. In 
order to define QCD beyond perturbation theory, Wil- 
son has regularized the theory on a 4-dimensional space- 
time lattice [H]. The quark fields then reside on the 
lattice points, while the gluon fields are naturally asso- 
ciated with the links connecting neighboring points. In 
lattice QCD the gluon field is represented by SU(3) ma- 
trices, which act as parallel transporters. Due to its non- 
perturbative nature and great complexity, lattice QCD 
requires large-scale Monte Carlo simulations. In this way, 
during the past decade very impressive progress has been 
made in determining the hadron spectrum as well as nu- 
merous static properties of hadrons. Also the nature of 
the transition that separates the low-temperature chirally 



broken phase from the quark-gluon plasma at high tem- 
peratures has been understood in some detail. However, 
despite of these successes, there are also important prob- 
lems in QCD that cannot be addressed with Monte Carlo 
methods. Such problems include QCD at large baryon 
density, which is relevant in nuclear and astrophysics, 
in particular, for understanding the state of matter in- 
side neutron stars, as well as all problems concerning the 
real-time dynamics of strongly interacting matter, in par- 
ticular, heavy-ion collisions. In these cases, severe sign 
problems prevent importance sampling which renders the 
Monte Carlo method inapplicable. 

Similar problems arise in strongly correlated systems 
in condensed matter physics. In recent years, several 
such systems have been emulated with ultra-cold mat- 
ter in optical lattices, thus realizing an analog quantum 
simulator. Since these systems realize quantum dynam- 
ics in their hardware, they do not suffer from the sign 
problem, and are thus ideally suited for investigating 
real-time dynamics. Recently, quantum simulator con- 
structions have been proposed for Abelian U(l) gauge 
theories with and without fermionic matter. Some of 
these constructions take advantage of the quantum link 
formulation of lattice gauge theory, which has a finite- 
dimensional Hilbert space per link. In Wilson's formu- 
lation, on the other hand, the Hilbert space is infinite- 
dimensional, and thus much harder to realize in ultra- 
cold matter. A non-Abelian quantum link model with 
gauge group SU(2) has been formulated in [SHE]- Sub- 
sequently, in [20l El] quantum link models have been 
constructed for the non-Abelian gauge groups SU(N), 
U(N), SO(N), and Sp(N), as well as for the exceptional 
group G(2). In particular, lattice QCD has been formu- 
lated as an SU(3) quantum link model [SO]- In order to 
reach the continuum limit, the model is formulated with 
an additional dimension of finite extent, which ultimately 
disappears via dimensional reduction. In this framework, 
quarks with a chiral symmetry emerge naturally as do- 
main wall fermions. While it may take a long time until 
a quantum simulator will be realized for quantum link 
QCD in the continuum limit, it is timely to investigate 
simpler non-Abelian gauge theories in order to address 
non-perturbative questions that are inaccessible to clas- 
sical simulation. 

In this paper, we construct quantum simulators for 
U(N) and SU(N) lattice gauge theories in (1 + 1) 
and (2 + 1) space-time dimensions coupled to staggered 
fermions with a Z(2) chiral symmetry. Being discrete, 
this chiral symmetry can break spontaneously in just 
one spatial dimension, at least at zero temperature. In 
two spatial dimensions it can break spontaneously even 
at finite temperature. In (2 + 1) dimensions, staggered 
fermions have an additional discrete Z(2) p flavor symme- 
try, which allows us to distinguish "mesons" from "glue- 
balls" . In a U (N) gauge theory, baryon number is part of 
the local symmetry. Hence, baryons violate Gauss' law 
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and thus do not belong to the physical spectrum. In an 
SU(N) gauge theory, on the other hand, there is a global 
£7(1) b symmetry that distinguishes physical states with 
different baryon numbers. In an SU(2) gauge theory two 
quarks arc confined inside a baryon. Then, in contrast to 
the real world, baryons are bosons. Despite numerous dif- 
ferences with QCD, concerning the space-time dimension, 
the gauge group, and the nature of the chiral symmetry, 
U{2) and SU{2) lattice gauge theories in (1 + 1) and 
(2 + 1) space-time dimensions are still useful toy models 
for investigating interesting non-perturbative dynamics. 
For example, one can study the real-time dynamics of 
chiral symmetry breaking across a phase transition, the 
superfluidity of baryons resulting from the spontaneous 
breakdown of U(1)b, and one can even mimic heavy ion 
collisions. 

A distinguishing feature of quantum link models is that 
the gauge fields can be realized as fermion bilincars. The 
fermionic constituents of gauge fields are known as ris- 
hons. As we will see, the rishon formulation is crucial for 
our construction of a quantum simulator for non-Abelian 
gauge theories. In fact, both "quarks" and rishons are 
represented by the same fermionic atoms that hop be- 
tween the sites of an optical lattice. From a theoretical 
perspective, the color indices of "quarks" and rishons can 
be contracted to form color-neutral "meson" , "glueball" , 
and "constituent quark" operators. It is remarkable that 
the theory can be formulated entirely in terms of these 
U (N) or SU{N) gauge invariant objects. Thanks to the 
explicit reduction to the gauge invariant sector, exact di- 
agonalization studies are possible even on moderate size 
systems, which can be used to validate an experimentally 
realized quantum simulator. Since the number of rishons 
per link is conserved locally, in the U (N) or SU (N) gauge 
invariant formulation of the theory, an additional U(l) 
link gauge invariance emerges. This ensures that the con- 
fining dynamics of "quarks" and "gluons" is still correctly 
represented by "mesons", "glueballs", and "constituent 
quarks" . 




TrUp 



H-O**- 



-Q « » l 



det U, 



FIG. 4. Rishon dynamics: The trace part of the Hamilto- 
nian induces a hopping of rishons of various colors around a 
plaquette. The determinant part shifts a color-neutral rishon 
baryon from one end of a link to the other. 



U(N) and SU(N) Quantum Link Models 

Let us consider a system of staggerd fermions on a d- 
dimensional spatial lattice coupled to a U(N) gauge field. 
The corresponding Hamiltonian then takes the form 

H = -t} j (S X y^ X U X y1py + ll.C.) + TJl Sj^f/'J. 

(xy) x 

2 /2 

i 9_ \ Ira t a , na r>a \ , 9 \ p,2 

2 ^ V-'xy^xy t ^xy^xy) ' ^ 2-^1 X V 

{xy) {xy) 

~ 4^2 X! {^U wx U X yU yz U zw + h.c.) 

(wxyz) 

- 7 (^tU xy + h.c.) . (7) 

(xy) 

Here and i\) l x (i € {1,2, ...,N} are fermion cre- 
ation and annihilation operators that obey standard anti- 
commutation relations, i is a hopping parameter, and 
Sx, Sxy are sign- factors associated with the points x and 
with the links connecting neighboring lattice points x 
and y, respectively. The site factor is given by s x = 
( — l) Xl ^ VXd . For the links in the 1-direction s xy = 1, 
for those in the 2-direction s xy = ( — l) Xl , and for the 
links in the fc-direction s xy = ( — l) TlH l- x *-i. JJ xy is 
an N x N matrix (with matrix elements Ul? y ) associ- 
ated with the oriented link connecting the neighboring 
points x and y, E xy , L% y , and R xy are Abelian and non- 
Abelian electric flux operators associated with the same 
link, and g' and g are the corresponding Abelian and 
non-Abelian gauge couplings. The magnetic field energy 
is represented by the term associated with the elemen- 
tary plaquettes (wxyz). Finally, the term proportional 
to 7 explicitly breaks a U (N) gauge symmetry down to 
SU(N). 

The various operators obey the commutation relations 

[L a ,L b ] - 2if abc L c , [R a ,R b ] = 2tf abc R c , 

[L a ,R b ] = [E, L a ] = [E, R a ] =0, 

[L a , U] = -\ a U, [R a , U] = UX a , [E, U] = U. (8) 

Here we have suppressed the link index xy. Operators 
associated with different links commute with each other. 
The Hermitean generators of SU(N) obey the commuta- 
tion relations 



[X a ,X b ] =2if abc \ c , Tv\ a \ b = 25 ab 



(9) 



where f a bc are the structure constants of the SU(N) alge- 
bra. By construction, the Hamiltonian of eq. Q is gauge 
invariant, i.e. it commutes with the infinitesimal genera- 
tors of SU (N) gauge transformations 



[G x ,G y ] — 2iS xy f abc G x 



(10) 



where A; is a unit-vector in the fc-direction. A general 
SU(N) gauge transformation is represented by the uni- 
tary transformation V = Y\ x exp(ia x G x ), which acts as 

U' xy = V^U xy V = exp«A a )*7„, exp(-ia^\ a ). (11) 

For 7 = the Hamiltonian has an additional U(l) gauge 
symmetry which is generated by 



[G x , Gx\ ~ 0. 



(12) 



A general U(l) gauge transformation, represented by the 
unitary transformation W — Y[ x Gxp(ia x G x ), acts as 

U'xy = W f U xy W = exp(ia x )U xy exp(-ia y ). (13) 

In Wilson's formulation of lattice gauge theory, the 
link matrices U xy are parallel transporters that take val- 
ues in the gauge group U(N) or SU(N), and E xy , L xy , 
and R xy are the corresponding canonically conjugate mo- 
mentum operators, which take derivatives with respect 
to the matrix elements of U xy . In that case, the com- 
mutation relations of eq.([8]) are realized in an infinite- 
dimensional Hilbert space per link, which poses great 
challenges for the implementation in a quantum simula- 
tor. Quantum link models extend the theoretical frame- 
work of lattice gauge theories by realizing the algebra of 
cq. ^ in a finite-dimensional Hilbert space per link. This 
is achieved by giving up the classical commutative nature 
of the matrix elements of U xy . In contrast to Wilson's 
lattice gauge theory, in quantum link models the real and 
imaginary parts of U x J y are no longer represented by real 
numbers, but by non-commuting Hermitean operators. 
In this sense, quantum links are gauge covariant general- 
izations of quantum spins. 

On each link the commutation relations of eq.Q can 
be realized by using the generators of an SU(2N) algebra. 
In particular, the real and imaginary parts of the N 2 ma- 
trix elements U x J y are represented by 2N 2 Hermitean gen- 
erators of the embedding algebra SU(2N). The 2(N 2 -l) 
generators L a and R a of SU(N) l X SU (N) r gauge trans- 
formations on the left and on the right end of the link also 
belong to the SU(2N) algebra. The Abelian 17(1) gauge 
transformations do not distinguish between left and right 
and are represented by the additional generator E. Al- 
together there are thus 2N 2 + 2(N 2 - 1) + 1 = AN 2 - 1 
generators, which form the algebra SU(2N). This re- 
quires that, unlike in Wilson's lattice gauge theory, the 
elements of a link matrix, U lJ , are non-commuting oper- 
ators that obey the commutation relations 

p/tf tfMt] = 2 (6 ik a%R a - 5 jlf jt k L a + 25 lk S 3l E), 

[U&Uu] = [(U f ) ij ,(U%}=0. (14) 

Again, we should emphasize that the commutation rela- 
tions are local, namely all commutators between opera- 
tors assigned to different links are zero. 



Since all elements of U xy commute with each other 
(although U xy and UJ. y do not commute) the definition 
of detU xy does not suffer from operator ordering ambi- 
guities. By construction, for 7^0, the Hamiltonian 
of eq.([7]) is invariant only under SU(N) but not under 
additional U(l) gauge transformations. However, there 
is a subtlety that needs to be taken into account. In 
the fundamental representation of SU(2N), the opera- 
tor that represents detU xy turns out to vanish. Hence, 
in that case even the Hamiltonian of eq.Q has a U(N) 
gauge invariance. On the other hand, since the above 
commutation relations can be realized with any repre- 
sentation of SU(2N), we can use higher representations 
in which the determinant in general does not vanish. 
The (27V)!/(A^!) 2 -dimensional representation of SU(2N) 
is the smallest with a non- vanishing determinant. For 
an SU(2) gauge theory the Hilbert space of the corre- 
sponding quantum link model then has 6 states per link, 
while for SU(3) it has 20 states per link. In the standard 
Wilson formulation of lattice gauge theory, on the other 
hand, the single link Hilbert space is already infinite- 
dimensional. 



Global Symmetries of the U(N) and SU(N) Quantum 
Link Model 



Let us discuss the global symmetries of the Hamilto- 
nian of eq.Q. First, we consider the shift symmetry Sk 
by one lattice spacing in the fc-direction 



(_-\ \Xk+i-\ t-Xd.l.i S h Tjij jjij 

K ' ^x+k 1 x v x+k.y+k 1 



Sk T a j a Sk T>a 

xy x+k,y+k' x y 



x+k^y+k 1 



Sk pa 
xy 



E a 

x+k,y+k ' 



(15) 



Only the mass term breaks this symmetry explicitly. 
Shifts by an even number of lattice spacings in the fc- 
direction play the role of ordinary translations. The as- 
sociated conserved quantum number is the momentum 
Pk €] — 5r, 5^]. In the zero- momentum sector the differ- 
ent shifts Sk form the group Z(2) d . Simultaneous shifts 
by one lattice spacing in an even number of different di- 
rections play the role of flavor transformations, which 
form the subgroup Z(2) d_1 . Simultaneous shifts by one 
lattice spacing in an odd number of different directions 
can be decomposed into a flavor transformation and a 
Z(2) chiral transformation, which is explicitly broken by 
the staggered fermion mass term. In one spatial dimen- 
sion, the chiral transformation is simply the shift by one 
lattice spacing. In two dimensions, simultaneous shifts 
by one lattice spacing in both directions represent a Z(2) 
flavor transformation, while a shift in either the 1- or the 
2-direction correspond to the Z(2) chiral symmetry. 

In a U(N) gauge theory, there are no baryons because 
the U(l) symmetry contained in U(N) is gauged. Hence, 
baryons (which are "charged" under this U(l) symmetry) 
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violate the Gauss law and are thus excluded from the 
physical Hilbert space. In an SU(N) gauge theory, on 
the other hand, there is a global U(1)b baryon number 
symmetry, i.e. [H, B] = 0, which is generated by 



N 



(16) 



Note that baryon number is counted relative to a filled 
Dirac sea. 

Another important symmetry is charge conjugation, 
which for staggered fermions also involves a shift by one 
lattice spacing. As for the chiral and flavor symmetries, 
there are separate symmetries Ck for the different direc- 
tions. 

Vx 1 > V x+k' U *V U x+k,y+k> 



C k j a 



-L/ a * 

x+k ,y+k 5 



h xy 



-R a * 

x+k^y+k ' 



-E a 

x-\-k,y+k ' 



(17) 

A combination of the two symmetries Ck and C'i acts as 

CiCk Tpi = - SiSk ^x- (18) 

This implies that, up to an overall minus-sign, the combi- 
nation of C'k and Ci corresponds to a flavor transforma- 
tion. A combination of an odd number of Ck operations, 
on the other hand, plays the role of a genuine charge 
conjugation. It should be noted that, strictly speaking, 
charge conjugation is not a Z(2) symmetry, because ap- 
plying it twice is equivalent to a translation, up to a 
minus-sign. Charge conjugation is a symmetry only if 
one chooses a non-chiral (i.e. real or pseudo-real) repre- 
sentation of the embedding algebra SU(2N). Finally, let 
us also list the parity symmetry P, which acts as 



V4 = 

T" 

L xy 



E X y — 



lb* P 

= R a _ 



w xy 



V,-xi 



U: 

^xy 



-y,-xi 



TO. 

L -y, 



-E- 



(19) 



are characterized by a lattice point x and a link direction 
±k. They obey canonical anti-commutation relations 



{ c x,±ki c il±i} - S xy 5±k,±i5ij, 

{^x,±k^ y ,±l} = { C x,±ki^y,±.li = °- 



(20) 



It should be noted that quark and rishon operators also 
anti-commute with each other. Under SU(N) gauge 
transformations the operators c and transform in the 
fundamental and anti-fundamental representation, re- 
spectively. It is straightforward to show that the com- 
mutation relations of eq.(|8| are satisfied when we write 

ra _ Jt \a J no _ if \a J 
^xy ^ij+^ij^i.+ i XL xy y,— ij y,— > 

1 



E r 



I 'A i 

2 \ y - v, 



c x ,+ c x,+)i Uj, — cj. i c£ . (21) 



All operators introduced so far (including the Hamilto- 
nian) commute with the rishon number operator 



Kf = c't J 



u x, + u x,+ 



(22) 



on each individual link. Hence, we can limit ourselves 
to superselection sectors of fixed rishon number for each 
link. This is equivalent to working in a given irreducible 
representation of SU(2N). Let us use the rishon repre- 
sentation to take a closer look at the determinant oper- 
ator that we used to break the U(N) gauge symmetry 
down to SU(N). We have 

Hpt/7 — — e ■ ■ TT i l i ' i n i2i ' 2 jf Nl ' N f-, ■/ •/ 

LLt^hUxy — AT .t llt2 ...z N U X y <~> X y ■ ■ ■ <J X y %>...%' 



1 

^%xi2—iN ^x, + u y, — u £c, + °y 
It „2 2t 



' • L a; : + L 2/,- Cl l l 2 



C x, + C y, — C x, + C y,- 



c N c N * 



(23) 



Only when this operator acts on a state with exactly 
ftf = N rishons (all of a different color), it can give a 
non-zero contribution. In all other cases the determinant 
vanishes. This means that we can reduce the symmetry 
from U(N) to SU (N) via the determinant only when we 
work with exactly Af = N fermionic rishons on each link. 
The number of fermion states per link then is 



The Rishon Representation of Quantum Links 

In this section we formulate quantum link models us- 
ing anti-commuting operators describing fermionic con- 
stituents of the gluons. The algebraic structure of a 
quantum link model is determined by the commutation 
relations derived above. The Hilbert space is a direct 
product of representations of SU (2N) on each link, with 
the generators on different links commuting with each 
other. We can thus limit ourselves to a single link. The 
commutation relations can be realized in a representa- 
tion using anti-commuting rishon operators c l x ±k , ±k 
with color index i £ {1, 2, N}. The rishon operators 
are associated with the left and right ends of a link and 



2N 
N 



(2N)l 
t/VTp' 



(24) 



This is the dimension of the SU(2N) representation with 
a totally antisymmetric Young tableau with N boxes (ar- 
ranged in a single column). 

The various symmetries act on the rishons as 



s k J 



-x,±l 

C k J _ 
c x,-l — 



C - k 

x+k,±V 



x,+l 



*t 



(25) 



- x+k -V " X > ±1 
It is straightforward to show that these transformations 



induce the corresponding transformations or eqs.(15l, 



(17), and (191 
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Glueballs, Mesons, and Constituent Quarks 

In this section we express the Hamiltonian of the U(N) 
quantum link model with quarks in terms of color neu- 
tral operators. Gauge invariance requires that these op- 
erators are local bilinear combinations of rishons and 
quarks, which we refer to as glueballs, mesons and con- 
stituent quarks. All these objects — including the con- 
stituent quarks — are bosons. The determinant term in 
the SU(N) quantum link Hamiltonian would give rise to 
an additional color-singlet rishon-baryon consisting of N 
rishons. For odd N this object would hence be a fermion. 
To avoid complications related to these objects we limit 
ourselves to U(N) in this section. 

By contracting the color indices of two rishons, we con- 
struct "glueball" operators 



x,±k,±l 



-x,±k c x,±h 



(26) 



which satisfy the local commutation relations of U (2d) 

[$x,±k,±h ®y,±m,±n] = 

= S xy (5±i,± m $ X: ±k,±n — 5±k,± n $x,±m,±l)- (27) 

By contracting the color index of two quark fields, we 
construct a "meson" operator 



m x = vM, 



(28) 



which generates a U(l) algebra. While glueballs and 
mesons are not related via their commutation relations, 
i.e. [& x ,±k,±i, My) — 0, they are both related with the 
constituent quark operators that one obtains by contract- 
ing the color indices of a quark and a rishon 



Qx,±k = cJ ±fc Vii 
via the commutation relations 



(29) 



[M x ,Qy t ± k ] = ~5 xy Q x . ±k - (30) 

Finally, the commutation relations of the constituent 
quark operators take the form 



[Qx,±kiQ y ,±i\ — $xy(&x,±k,±l — $±k,±lM x ), 



[Qx,±k,Qy,±i] - [Q x ,±k>Qv,±i] = °- 



it 



(31) 



Thus, the inclusion of the constituent quark operators 
completes the site-based algebra of U(2d) x U(l) to 
U(2d+l). 

Let us now consider a single rishon per link, i.e. Af xy — 
1. Then detU vanishes and the electric flux terms in 
the Hamiltonian of eq.([7| (proportional to g 2 and g' 2 are 
trivial constants. Expressing the Hamiltonian in terms 
of glueball, meson, and constituent quark operators then 



leads to 



(my) 



I'/ 



2 ^ . (UwxUxyUy Z U zw -\~ ll.C.) 



(wxyz) 

tJ2 ( s xyQi t+k Qy,-k + h.C. 



y^^s x M x 



(xy) 



I'/ 



+l,+k < &y-k,+l < bz,-l,-kl32) 



{wxyz) 



where y — x + k, z = y + I, and w = x + I. Remark- 
ably, this expresses the Hamiltonian entirely in terms 
of U(N) gauge invariant operators. Still, it should be 
noted that the number of rishons per link is conserved, 
i.e. [H,Af xy ] = 0. Also note that 



(33) 



This framework is well suited for exact diagonalization 
studies because it explicitly eliminates gauge variant 
states. 

Under the various global symmetries the glueball, me- 
son, and constituent quark operators transform as 



®x,±l,±m — ® x+ k ,±l,±m' 

M_ l£ , Sk Q x ,±i = Q 



M 



x+k , ^ x ,±t ^x+k^V 
> ' '" = — S ± tS ± m (^x+k,±m,±l 



NS 



±l,±m) 



Ch M x 

Ck Q x ,±i = T (-i) Xl+ - +Xk Q f 



x+k,±l' 



?x,±l,±m — ^-i^i^m, 
Pus _ n.r Pi 



M x — M- x , Q x ,±i — Q-xrfi 



(34) 



Here the sign s±i = ±1 is the sign in ±Z. The above 
symmetry transformations follow directly from the trans- 
formation properties of the fermions and of the rishons 



(eq.(25l). It is straightforward to show that they leave 



the commutation relations between the glueball, meson, 
and constituent quark operators invariant. 



SUPPLEMENTARY MATERIAL II: 
NON-ABELIAN LATTICE GAUGE THEORIES IN 
ALKALINE-EARTH GASES: THE EFFECT OF 
IMPERFECTIONS IN ATOMIC 
IMPLEMENTATIONS 

Introduction 

Imperfections in atomic physics implementations have 
to be carefully analyzed in order to understand the mini- 
mal experimental requirement to perform accurate quan- 
tum simulations. Here, we discuss how typical experi- 
mental errors affect the precision of quantum simulations 
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of non-Abelian lattice gauge theories in alkaline-earth se- 
tups, considering the paradigmatic case study of U{2) 
quantum link models. In the next section, we briefly 
describe our model Hamiltonian, its experimentally rel- 
evant perturbations, and a set of observables quantita- 
tively describing the departure from gauge invariance. 
In the last section, we present numerical results based on 
Lanczos diagonalization of small systems to underpin the 
validation issue at a quantitative level. 

Microscopic model Hamiltonian, imperfections, and 
observables 



in the U (2) context. In the following, we particularly fo- 
cus on the effect of different tunneling rates for the two 
color species, investigating to what extent they can affect 
gauge invariance. 

Accuracy of gauge invariance: relevant observables 

In this section we investigate how deviations from a 
fixed rishon number per link in the microscopic Hamil- 
tonian H and color-dependent tunneling rates in -ffi mp 
affect the accuracy of the gauge symmetries. 



For simplicity here we consider a (1 + 1)D U(2) QLM 
with two rishons per link. As discussed in the main text, 
the QLM Hamiltonian H is induced in second order per- 
turbation theory from the microscopic Hamiltonian 

H = UJ2(^+ + n x+ i,- - 2) 2 + m^(-l)^^i 

X X 

+ tJ2(c x \+i>l + 4 f ,_C + h.c), (35) 

X 

where n x ± = c x± c x± . This Hamiltonian is manifestly 
SU(2) gauge invariant, i.e. [H, G%] = 0. Unlike the QLM 
without the detU term, it is, however, not necessarily 
also U{\) gauge invariant. Also the constraint of J\f xy = 2 
rishons per link will only be satisfied approximately for 
sufficiently large U. 

In an experimental realization with ultracold alkaline- 
earth atoms in an optical lattice, one may expect 
some deviations from the ideal microscopic Hamiltonian. 
While interactions do not need any particular fine-tuning 
in our implementation, since the difference in scattering 
lengths for the different mj states are negligible com- 
pared to all other energy scales, different tunneling rates 
and local potentials require accurate lattice design, and 
may be subject to experimentally relevant imperfections. 
Such imperfections may give rise to color-dependent tun- 
neling rates £j. The actual Hamiltonian of a slightly im- 
perfect experimental realization may take the form 

H imp = Uj2(^++n x+1 ^-2) 2 + mJ2(-^ X ^^ 

X X 

+ £t i (< + ^ + 4_^+h.c.) 

X 

+ 2AC/^(n X;+ -l)(n x+1 ,_-l). (36) 

X 

It should be noted that repeated indices are still summed. 
Due to the color-dependent hopping terms U, this Hamil- 
tonian is no longer SU(2) gauge invariant. The AU term, 
resulting from different overlaps of the Wannier func- 
tions in the partially overlapping rishon-sites, provides 
a source for the implementation of the detU term in the 
SU(2) model, but should be viewed as an imperfection 



Abelian U(l) gauge symmetry 
Let us first discuss the U(l) symmetry generated by 
G x = nl+nl - ^(K,x+i+Mx-i,x), (37) 

where 

ni=i/>M + + $,-<- ( 38 ) 

counts the number of fermions of color i on a basic build- 
ing block consisting of a quark- and the two adjacent 
rishon-sites forming a triple-well. Since the fermion num- 
bers n x are exactly conserved by construction, one should 
check whether the number Af x ,x+i = 2 of rishons per link 
remains constant. Deviations from this constraint are 
quantified by 

(Afx.x+i) =2 + Aiink, (39) 
which should be small for sufficiently large U. 

Non-Abelian SU(2) symmetry 

We now focus on the non-Abelian part of the symmetry 
group, SU(2), and check the expectation values of the 
generators 

G a x = ^otM + 4+^4,+ + (40) 

which should be zero in case of exact gauge symmetry. 
Here, a a are the Pauli matrices. Notice that the diagonal 
U(l) subgroup of the SU(2) symmetry is also automati- 
cally fulfilled. This is because the generator 

Gl=nl- r? x (41) 

always has a vanishing expectation value by initial state 
preparation: the number of particles assigned to each 
basic building block (consisting of a quark- and two ad- 
jacent rishon-sites forming a triple-well) is strictly con- 
served. Here, we focus on the generator G\ and calculate 
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FIG. 5. Deviations from the conserved number of rishons as 
a function of t/2U for different values of m/t. 



the expectation value in the ground state of H lmp (with 
AU = 0) 

A^E^) 2 )- ( 42 ) 

X 



Numerical results 

In this section, we present numerical results based on 
the previous theoretical considerations. 

U(l) rishon number symmetry: validation 

In Fig. [5] we present exact diagonalization results for 
a system with 12 sites (4 quark- and 8 rishon-sites) . This 



shows how the number of rishons per link, indicating the 
accuracy of the U(l) symmetry, is very close to the de- 
sired value even for relatively small values of the interac- 
tion strength U . This is similar to the previously studied 
U(l) case. 
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TABLE II. Discrepancy from exact gauge symmetry in the 
microscopic model as a function of z for different inter-particle 
interactions U. 



SU(2) symmetry: validation 

In Table [TTJ we list typical values of D\ as a func- 
tion of the interaction strength U (still considered SU(2) 
symmetric) for approximate tunneling symmetries, where 
£2 =t,ti = (1— z)t. Remarkably, even for relatively large 
values of z ~ 0.1, errors due to approximate gauge invari- 
ance are extremely small, well beyond all other possible 
imperfections such as, e.g., the effect of the trapping po- 
tential or the space-dependent tuning of optical Feshbach 
resonances. 



